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Abstract 

In this paper, we investigate the linear perturbations of the spheri- 
cally symmetric spacetimes with kinematic self-similarity of the second 
kind. The massless scalar field equations are solved which yield the 
background and an exact solutions for the perturbed equations. We 
discuss the boundary conditions of the resulting perturbed solutions. 
The possible perturbation modes turn out to be stable as well as un- 
stable. The analysis leads to the conclusion that there does not exist 
any critical solution. 
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1 Introduction 

In General Relativity, gravitational collapse of a realistic body is one of the 
important problems. At the threshold of black hole formation, the matter is 
just about to form a black hole and an infinitesimally small perturbation can 
either cause the matter to disperse to infinity or to form a black hole. The 
dynamics close to the threshold exhibits interesting behavior such as power 
law scaling of length scales, self-similarity of the solutions and universality. A 
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typical setup is to take an initial matter distribution parameterized by a single 
parameter. The critical solution is self-similar, i.e., it repeats itself on ever 
decreasing length scales. The fluid density and pressure increase by many 
order of magnitudes and in some cases, the fluid velocity becomes extremely 
relativistic. Another way of investigating the critical collapse is to use the 
fact that the critical solution is continuously self-similar. Self-similarity leads 
to the set of partial differential equations into a set of ordinary differential 
equations that can then be solved quite easily numerically to a very high 
precision. One then performs a perturbation analysis around the critical 
solution to determine the scaling exponent. 

Critical phenomena in gravitational collapse were discovered by Choptuik 
[1,2] in the spherically symmetric collapse of massless scalar field. Since then 
the phenomenon has been observed in a variety of matter sources. There are 
two steps to find critical solutions: Firstly, find a generic family (or families) 
of solutions, defined by a parameter, say p such that when p > p* the collapse 
forms black holes and when p < p*, it does not. Secondly, the perturbations 
of the solution p = p* are performed to investigate the spectrum of their 
modes. If the solution has only one unstable mode, then this solution is a 
critical solution and the exponent 7 is given by 




where a is the unstable mode [3]. 

Garfinkle [4] found a class, say, S[n], of exact solutions to the Einstein 
massless scalar field equations in (2 + l)-dimensions. He showed that in the 
strong field regime the n = 4 solution fits very well with the numerical crit- 
ical solution found by Pretorious and Choptuik [5]. Later, Garfinkle and 
Gundlach [6] studied their linear perturbations and found that only the so- 
lution with n = 2 has one unstable mode, while the one with n = 4 has 
three. They further required that no matter field should come out of the 
already formed black holes. This additional condition seems physically quite 
reasonable and has been widely used in the investigation of black hole pertur- 
bations [7]. Hirschmann et al. [8] systematically studied the critical gravita- 
tional collapse of a scalar field. They surveyed all the analytic, continuously 
self-similar solutions and also examined their perturbations considering their 
global structure. Clement and Fabbri [9] investigated analytical treatment 
of critical collapse in 2 + 1-dimensional AdS spacetime. Cavaglia et.al. [10] 
analysed approximately self-similar critical collapse in 2 + 1-dimensions. 
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Miguelote et al. [11] studied the gravitational collapse of self-similar per- 
fect fluid in 2 + 1-dimensional spacetimes with circular symmetry. They also 
studied the linear perturbations of homothetic self-similar stiff fluid solutions. 
Frolov [12] studied the perturbations of the continuously self-similar critical 
solution of the gravitational collapse of a massless scalar field. Brandt et 
al. [13] studied the gravitational collapse of spherically symmetric perfect 
fluid with kinematic self-similarity. They also studied the solutions of the 
Einstein field equations found by Benoit and Coley [14] and concluded that 
some of the solutions represent gravitational collapse. Wang [15] studied 
the critical collapse of a cylindrical symmetric scalar field. He introduced 
the notion of homothetic self-similarity to four-dimensional spacetimes and 
then presented a class of exact solutions to the Einstein massless scalar field 
equations. Wang et al. [16] studied plane symmetric self-similar solutions 
to Einstein's four-dimensional theory of gravity and explored the local and 
gravitational conditions. 

Recently, Chan et al. [17] investigated the solution of the Einstein mass- 
less scalar field equations with kinematic self-similarity of the second kind in 
the (2+l)-dimensional spacetimes with circular symmetry. They discussed 
their local and global properties and also found that some of these solutions 
represent gravitational collapse of the scalar field. The same authors [18] 
studied the linear perturbations of the (2+l)-dimensional circularly symmet- 
ric solution with kinematic self-similarity of the second kind. They obtained 
an exact solution for the perturbation equations and the possible perturba- 
tion modes and showed that the background solution is a stable. 

In this paper, we investigate the linear perturbations of the spherically 
symmetric spacetimes with kinematic self-similarity of the second kind. We 
analyze these solutions to see whether the solution is stable or not. The pos- 
sible perturbation modes are also discussed. The paper has been organized as 
follows. In the next section, we shall consider the general spherically symmet- 
ric spacetimes and calculate self-similar variable of the second kind. Section 

3 is devoted to the linear perturbation of the field equations. In sections 

4 and 5, we find the possible solutions of the linear perturbation equations 
and check the boundary conditions for the perturbed solutions respectively. 
Section 6 is focussed for discussion of the results. 
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2 Self-Similarity of the Second Kind 



In this section, we calculate the self-similar variable of the second kind for 
spherically symmetric spacetimes. Further, we write down the non-zero com- 
ponents of the Ricci tensor in terms of the self-similar variable. The most 
general form of the spherically symmetric metric is given by [19] 

ds 2 = e 2 ^dt 2 - e m ^dr 2 - r 2 S 2 (t,r)(d6 2 + sin 2 M/> 2 ), (2) 

where $(t,r), fy(t,r), S(t,r) are arbitrary functions of t, r. The non-zero 
components of the Ricci tensor are 

i? 00 = e 2 ^[$ r ($ r -M/ r + ^ + -) + $ rr ] 

b r 

- ^ + ^ + $^-* 2 -^, (3) 
2% 2^ t S r 2<Z> r S t 2S t 2S rt 

Roi — 1 c 1 c c c - ' ^) 

r b b rb b 

R u = e 2(*-*)[^ + ^(^ t + _££ -$ t )]-$ rr 

^ 2 2S rr 4S r 2^ r S r 2V r 
b rb b r 
R22 = e- 2 *(r 2 S? + r 2 SS tt -r 2 SS t <l> t + r 2 SS t y t ) 

+ e-^i-S 2 - r 2 S 2 - r 2 SS rr + rS 2 V r + r 2 SS r V r 

- 4rSS r -rS 2 <5> r -r 2 SS r <5> r ) + l, (6) 
R 33 = R 22 sm 2 6, (7) 

where the subscripts t and r mean differentiation w.r.t t and r respectively. 
We consider here solutions with kinematic self-similarity of the second kind 
for which the self-similar variable x turns out to be 

x = ln(— ^-r), T = -ln(-t), (8) 

{-t)<* 

or inversely 

ax — T 

r — e a , t — —e , (9) 

where a is a dimensionless constant. The components of the Ricci tensor in 
terms of the self-similar variable can be written as follows 

e 2(*-tt) g l 

R ™ = —^2— {^x-*x + 2-^ + 1) + <f> xx } -—{ a 2 [m TT 
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+ * T (1 + ^ r - $ T ) + 2^ + 2^(1 - $ T )] + a[2V TX + * T (* a 

- $ x ) + vj/ x (* r - $ r + 1) + 4^ + 2^(1 - $ T ) - 2^] 

+ [*« + *x(*,-$x)+2^-2^]}, (10) 

R« = -^d + |) + |(^-D-^] 

+ **(l + f) + f(*x-l)-^}, (11) 
2(*-*) C 

i?n = -^^{a 2 ^ rr + ^ T (l + ^ T -$ r + 2^)] 

+ a[2^ r:c + * X {1 + 2^ T - $ T + 2^) + * T (-$ X + 2^)] 

+ # xx + # x (* x - $ x + 2%)} + - ®* + 1) 

- $ xx + 2* x (l + ^)-|(S xx + S x )], (12) 

-2$ C 

R22 = r 2 S 2 {^^[a 2 (S T (l + ^ T -$ T + ^) + S TT ) 

+ a(S r (* x - $ x + 2^) + 5,(1 + ^ r - $ T ) + 2S r:c ) 

+ <? x (* x -<t> x + ^) + S xx ]- — [ -(S x (3 + $ x 

- *x + y) + S M ) + ^-* x + l]} + l, (13) 
#33 = i? 22 sin 2 6>. (14) 



3 Linear Perturbation of the Field Equations 

This section is devoted to set up the Einstein field equations by using linear 
perturbation. For this purpose, we take 

$(r,x) = $ (x) + e$x(a;)e fcT , 

yj/( T ,x) = V (x) + e^ l (x)e kT , 

S(t,x) = S (x) + eS 1 {x)e kT , 

(p(r,x) = cf) (x) + e<p l (x)e k \ (15) 
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where e is a very small real constant and k is an arbitrary constant. The 
quantities with subscripts and 1 denote background self-similar solutions 
and perturbations respectively. It is understood that there may be many 
perturbation modes for different values (possibly complex) of the constant 
k. The general perturbation will be the sum of these individual modes. The 
modes with Re(k) > grow as r — > oo and are referred to as unstable modes 
while the ones with Re(k) < decay and are referred to as stable modes. By 
definition, critical solutions will have one and only one unstable mode. 
We take the following background solution [18] 

$o(x) = 0, V (x) = -^ax, 

S (x) = 7^—e- l i ax , <f>o(x) = 2qhi(-t) (16) 
and the apparent horizon is given by 

r AH (t) = [(2-a)V=t}^, a<2, (17) 

where q = ±^g- For a massless scalar field 0, the Einstein field equations 
are 

R a b = K(j) a (t) b: a, 6 = 0,1,2. (18) 

Here we choose units such that k — 1 for the sake of simplicity. Using Eq.(15) 
in Eqs.(10)-(14), it follows that 

Rab = Rab(T,x,e) (19) 

If we take R a b as a function of e only and expand it in terms of e, it follows 
that 

Rab(r,x,e) = -J-^{R^( x ) + eR^(x)e^ + 0( e 2 )}. (20) 

The non-vanishing components of the Ricci tensor upto first order in e are 

R ( £(x) = e 2 ^°-*°- 3:+ «){$( ) (2$ / 1 + 

+ $; + 2($x - tfiMCS'o - %) + $o + & ] + ^[-2^<t> 

Oo oo 

+ 4($! - #i)$' S + 2%S[ + 2^S' ]} + ^{^(afc*! 
+ *;) + %(ak§ 1 + $;) - 2^' (aA;^ 1 + - a 2 k 2 V 1 
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- 2ak^[ - - a 2 k 2 ^! l - + —[2<l>' ( y akS 1 + S[) 

>->o 

+ 2S' (ak$ 1 + $i) - 2a 2 A; 2 S'i - 2a 2 fc^ - 4afcS; - 2^' 

- 2a 2 kS 1 - 2aS[ - ^(2S' <S> - 2S' ' - 2c*S£)]}, 

R$(x) = -2^^[-^(s' -^ s' -%s' + s' f ) 

ab d 

- $' (afc5i + S-i) - tt' S; - S^afcWi + ^ + $i) 
+ aifcSj + S'( - SoiahVi + ^) + akS 1 + 

RS(x) = 2 [2(^i-$i)« + ^ + ^o 

ar 

- $' * + 2 %^0) + 2«A;* , * 1 + 2#' *i + a 2 ^ 
+ 2afc*i + + a 2 k^ 1 + - a£;$ *i - ^' % 

- ak<f>!% - $'^0 + ^(akS S 1 % + 5 5^o 

<->o 

+ $'/ - $; + 2$' $' 1 - $ *; - J-(d^ + S'^ 

D 

+ W / 1 -^-^) + 2§(W-^-^ )], 

D 

+ 5 ' + 3S' + S + ^) + & S S[ - %S S[ + 2S' S[ 
D 

+ 3S S[ - ^S'oS, + ^S'^ - 35^1 + $^0^0 

- %S S' + S Q S'I - S^S, + <S>[S 2 " " 

Do 

+ e 2{ ^r T+x -* o) {ak(aS Q S 1 + akS S l + 2d 1 d\ - $' 5 5i 
+ d-o^i - $i^o + *iS' + 250) + «^o^ + 2S' S[ 

+ S S[% - S S[& + SqS'^ - SoS'^ + S S' ' - 

DO 

+ (5i - 2$ 1( S )(«5o + S' % - S' & + SZ + f^)}, 

D 

= R$(x)sm 2 9. 



Now we can calculate the quantities A ab = 4> a 4>b using Eq.(15) in Eqs 
(14). Thus 

A ab (r, x, e) = -fp^{4$(x) + eA^(x)e^ + 0(e 2 )}. 
The perturbed part is given by 



A {1 \ 

^00 1 






e 2r 






a 








A'S< 




= o, 








= 




^■33 1 




= 0. 





The linear perturbation equations can be written as 

BSt(x)=A2j>(x). 
Using Eqs. (27) and (28), it turns out that 

3 3 

Aaq{ak(f>i + </)[) = -a 2 k$i + -a^[ + a 2 k 2 ^ 1 + 2ak^[ 

+ V / ; + ?—^e 1 2 ax [(2k 2 + 2k + l)a 2 S 1 
+ 2a(2fc + l)S'i + 2^], 

o = 

2aq<j)' 1 = (a - 2)(ak^ 1 + %) + 

9 — n i 
+ -^e3 ra [a(2fe + l)Si 

+ (2ak + a + 2)S[ + 2S"\, 
= 2a 2 k(k-l)^ 1 + 2a(2k-l)^' 1 
+ 2^/ + a{ak^ + $[) + a 2 (^ 1 - $i) 

- a(2 ~ a) e^ 3! [a(2fc+ 1)^ + 2^], 

o = 2 (2 - «) + (« -2) + 

- (2-a)e5 a:E [ ai Si + (a + 2)^ + 2^ / ], 
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= a(uk<§ 1 + $[) - a(ak^ x + %) 

- a 2 $x -^^[a 2 ^ - l)5i 

+ a(4k-l)S[ + S'l], (34) 

= (a-2)(V 1 -& 1 )-V 1 (a 2 -4a 

2 — n 1 

+ A)-—e^ x [{A-a)S l 

+ {§-a)S' l + 2S'l\. (35) 



4 Solutions of the Linear Perturbation Equa- 
tions 

Now we solve system of the perturbed Eqs.(29)-(35). From Eq.(30), we have 

$i = b, (36) 

where b is an integration constant. Multiplying Eq.(35) by 2, adding in 
Eq.(33) and using Eq.(30), we obtain 

(a - 2)^i + 2e* ax [S 1 + 2S[ + S'l\ = 0. (37) 

Using Eq.(36) in Eq.(34), it becomes 

= a 2 (k-l)b-a%-a 2 k^ 1 
2 — rv i 

+ — — e2 Q > 2 A;(2fc-l)S 1 + a(4A;-l)^ + ^ / ]. (38) 
From Eq.(37), we have 

#i = ^— e* ax [S 1 + 2S' 1 + S'{\. (39) 
2 — a 

Eliminating and ty[ from Eqs.(38) and (39), it turns out that 

ASi + BS[ + CS'l + DS'l' + Ee"^ = 0, (40) 

where 

A = 2a 2 (2k+l)-a 2 k(2-af(2k-l), (41) 

B = 2a(4aA; + 2a + 2)-a(4A;-l)(2-a) 3 , (42) 

C = 2a(2ak + a + A) - (2 - a) 3 , (43) 

D = 4a, (44) 

E = -2ba 2 (k-l)(2-a). (45) 
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It can be found that the solution with self-similarity of the second kind is 
identical to the solution with first kind for the same type of fluid. This means 
that the spacetime can have two different kinds of self-similarities, i.e., there 
exist two vector fields where £^ describes self-similarity of the first 

kind and £^ of the second kind. This happens when the spacetimes has high 
symmetry. As the background solution with self-similarity of the second kind 
becomes identical to the solution with self-similarity of the first kind, we take 
a = 1. Thus Eqs.(41)-(45) become 

A = -2k 2 + 5k + 2, (46) 

B = 4k + 9, (47) 

C = 4k + 9, (48) 

D = 4, (49) 

E = -2bk + 2b (50) 

and Eq.(40) takes the form 

(-2k 2 + 5k+2)S 1 + {Ak + 9)S' l + {4k + 9)S'l + 4S';' + (-2bk+2b)e-^ x = 0. (51) 

This has the following solution 

where C\ is an integration constant and 

U = —P ^7 + —A, (53) 

12 12A 12 ' y J 

p = -9 -4k, (54) 
B = 27- 24k- lQk 2 , (55) 

A = [297- 756fc + 288fc 2 -64jfc 3 + 6v / 3v / A]3, (56) 
A = 999 - 4644A; + 6984A; 2 - 3936A; 3 + 1600fc 4 - 512A; 5 . (57) 



5 Boundary Conditions for the Perturbed So- 
lutions 

In this section, we shall discuss some geometrical and physical conditions 
[20] needed to be imposed for the spherically symmetry. For gravitational 
collapse, we impose the following conditions: 
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(i) There must exist a symmetry axis which can be expressed by 



X = ^\^ e) e ie) 9a h \ - (58) 

as r — > 0, where we have chosen the radial coordinate such that the axis is 
located at r = and £^ is a Killing vector with a close orbit which is given 
by ££ e) d a = d e . 

(ii) The spacetime near the symmetry axis is locally flat which can be written 
by [19] 

X a X h g ab - -1 (59) 

as r — > 0. It is mentioned here that solutions failing to satisfy this condition 
may also be acceptable. Since we are mainly interested in gravitational col- 
lapse, we assume that this condition strictly holds at the beginning of the 
collapse so that we can be sure that the singularity to be founded later on 
the axis is due to the collapse. 

(iii) Closed timelike curves can be easily introduced in spacetimes with spher- 
ical symmetry. For their absence the condition, 

Z(o$e)9<* < 0, (60) 

must hold in the whole spacetime. In addition to these conditions, we also 
require that the spacetime is asymptotically flat in the radial direction. For 
self-similar solutions, this condition cannot be satisfied unless we restrict 
their validity only up to a maximal radius say, r = r (t), and join them 
with others in the region r > ro(t) which are asymptotically flat as r — > oo. 
Here we simply assume that the self-similar solutions are valid in the whole 
spacetime. 

The boundary conditions at the event horizon tah, given by Eq.(17), take 
the following form 

r A H = ~t (61) 
and the corresponding metric becomes 

ds 2 AH = -4(-t) 2 (d6 2 + sin 2 0d<t> 2 ). (62) 

This shows that the apparent horizon is singular only at t — and the final 
state of the collapse is marginally naked singularity. Now we discuss the 
following three cases: 

(i) A > 0, (ii) A = 0, (iii) A < 0. 
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Case (i): In this case, using Eq.(52) in (39), we have 



For the first boundary condition, given by Eq.(58), we calculate the quantity 
\[X = rSi by using Eq.(8) in (52) so that 

. r . tj 8b(k — 1) , r , i , . 

S > = ^» - M»-lfflfc + 3 '^» (64) 
r u+1 8Wk-l) , 

rS ' = c 'HF" 8F-i 6fc ; 3 ( - ri) " < 65) 

For rSi — > as r — > 0, all the exponents of r must be greater than zero. It 
can be checked from Eqs.(53)-(56) that U + 1 is positive for — oo < k < 1 but 
turns out to be complex otherwise. Thus the first condition is satisfied for 
— oo < k < 1 but is not fulfilled otherwise. Consequently, these perturbations 
are limited by the boundary conditions for 1 < k < oo. 
Case (ii): For A = 0, k -> 1.70872 along with Eqs.(52) and (39), we have 

#i = 2.88766 + 0.2043e" a8196:E . (66) 

The boundary conditions, given by Eqs.(58)-(60), can be applied by making 
use of Eq.(8) in Eq.(52) 

S 1 = (-^)- L3196 + 5.77526(-^)-5, (67) 

rS 1 = r- a3196 (-t) L3196 + 5.77526r5(-t)i (68) 

Since all the exponents of r are not positive, the condition given by Eq.(58) 
is not satisfied. Again perturbations are limited by the boundary conditions. 
Case (iii): Here we use Eq.(52) in Eq.(39) so that 

11/1 = 8fc2 4 -l6fc + 3 + (1 + f/l)ci(coS ^ 2 + isinf/ 2)e^, (69) 



where 



= i(-9-«)- 27-24t-16f 



12 12(297 - 756/c + 288P-64A; 3 -108A)< 
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" ^4 arctail ( 297 _ 7 ' 6k I - 64k* )]} 

+ — (297- 756/c + 288/c 2 -64A; 3 -108A)s 



L 3 V 297- 756A; + 288A; 2 -64A; 3/J 

27 - 24k - lQk 2 

U 2 - {[ 



12(297 - 75Qk + 288F - 64P - 108A)s 
+ -U297 - 756A; + 288fc 2 - 64A; 3 - 108A)s] 



. .1 , 

x sin - arctanl — -r )\Yx, 71 

L 3 V 297- 756A; + 288P-64A; 3;JJ ' v ; 

r 1 , n r 27 - 24A; - 16A; 2 

^3 = { — (-3-4fc)-[ 



12 12(297- 75Qk + 288P - 64P - 108A) 5 

- -^-(297- 756A; + 288fc 2 -64A; 3 -108A)s] 

r i . ev^v/^A 

x cos — arctanl — — ) \\x. (72) 

L 3 V 297- 756A; + 288A; 2 -64A; 3;JJ v ; 



For the boundary condition, given by Eq.(58), we require that VX = rS\ 
and S\ is real only if c\ = 0. Thus we obtain 

51 " "8P-16A; + 3 ( ^ ) ' (73) 
86(fc-l) , ,1 

Clearly, rSi — > as r — > 0, the condition, given by Eq.(58), is satisfied only 
if k 7^ 1± ^/|. For the second boundary condition, given by Eq.(59), we have 

v v ab 81926 4 (/c - l) 4 r V 1 . 
Thus X a X b g ab -> as r -> 0. 
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6 Concluding Remarks 



We have investigated the linear perturbations of the spherically symmetric 
spacetimes with kinematic self-similarity of the second kind. The self-similar 
variable and the background solution found for the spherically symmetric 
spacetimes become identical to that of the circularly symmetric metric [18]. 
However, the linearly perturbed solution obtained for the spherically sym- 
metric metric is different from the exact solution for the circularly symmetric 
metric. The boundary conditions for all possible values of A are discussed 
for this linearly perturbed solution. 

For A > and A = 0, the first boundary condition, given by Eq.(58), is 
not satisfied. This shows that the perturbations are limited by the boundary 
conditions in both these cases. For A < 0, the boundary conditions are 
satisfied only for k ^ 1 ± which admits both stable and unstable modes 
for the perturbation. The stable modes are for k < and the unstable modes 
are for k > but k ^ 1 ± \J\. The unstable modes for the perturbation imply 
that it is not a critical solution. 

Acknowledgment 

The author would like to thank referee for the constructive comments. 

References 

[1] Choptuik, M.W.: Phys. Rev. Lett. 70(1993)9. 

[2] Choptuik, M.W., Hirchmann, E.W., Liebling, S.L. and Pretorius, S.: 
Phys. Rev. D68(2003)044007. 

[3] Hara, T., Koike, T. and Adachi, S.: Phys. Rev. D59(1999) 104008. 

[4] Garfinkle, D.: Phys. Rev. D63(2001)0044007. 

[5] Pretorius, F. and Choptuik, M.W.: Phys. Rev. D62(2000) 124012. 

[6] Garfinkle, D. and Gundlach, C: Phys. Rev. D66(2002)044015. 



14 



[7] Chandrasekhar, S.: The Mathematical Theory of Black Holes (Clarendon 
Press, Oxford University Press, Oxford, 1983). 

[8] Hirchmann, E.W., Wang, A. and Wu, Y.: Class. Quantum Grav. 
21(2004)1791. 

[9] Clement, G. and Fabbri, A.: Class. Quantum Grav. 18(2001)3665; 
Nucl. Phys. B630(2002)269. 



[10 
[11 

[12 
[13 

[14 
[15 
[16 
[17 

[18 

[19 

[20 



Cavaglia, M., Clement, G. and Fabbri, A.: Phys. Rev. D70(2004)044010. 
Miguelote, A.Y., Tomimura, N.A. and Wang, A.: Gen. Rel. Grav. 



36(2004)1883; arXiv: gr-qc/0505062 



Frolov, A.V.: Phys. Rev. D56(1997)6433; ibid 59(1999)104011. 

Brandt, C.F.C., Lin, L.M., Villas da Rocha, J.F. and Wang, A.: Int. J. 
Mod. Phys. Dll(2002)155. 

Benoit, P.M. and Coley, A.A.: Class. Quantum Grav. 15(1998)2397. 
Wang, A.: Phys. Rev. D68(2003)064006. 

Wang, A., Wu, Y. and Wu, Z.C.: Gen. Rel. Grav. 36(2004)1225. 

Chan, R., da Silva, M.F.A., Villas da Rocha, J.F. and Wang, A.: Int. 
J. Mod. Phys. D14(2005)1049. 

Chan, R., da Silva, M.F.A. and Villas da Rocha, J.F.: Int. J. Mod. 
Phys. D14(2005)1725. 

Stephani, H., Kramer, D., MacCallum, M., Hoenselaers, C. and Herlt, 
E.: Exact Solutions of Einstein Field Equations (Cambridge University 
Press, 2003). 

Carlip, S.: Quantum Gravity in (2+1) -Dimensions (Cambridge Uni- 
versity Press, Cambridge, 1998). 



15 



